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Abstract 



We present the second of two articles on the small volume fraction limit of a nonlocal 
Cahn-Hilliard functional introduced to model microphase separation of diblock copoly- 
mers. After having established the results for the sharp-interface version of the functional 
( [5] ) , we consider here the full diffuse-interface functional and address the limit in which e 
and the volume fraction tend to zero but the number of minority phases (called particles) 
remains 0(1). Using the language of T-convergence, we focus on two levels of this conver- 
gence, and derive first- and second-order effective energies, whose energy landscapes are 
simpler and more transparent. These limiting energies are only finite on weighted sums 
of delta functions, corresponding to the concentration of mass into 'point particles'. At 
the highest level, the effective energy is entirely local and contains information about the 
size of each particle but no information about their spatial distribution. At the next level 
we encounter a Coulomb-like interaction between the particles, which is responsible for 
the pattern formation. We present the results in three dimensions and comment on their 
two-dimensional analogues. 

Key words. Nonlocal Cahn-Hilliard problem, Gamma-convergence, small volume- 
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1 Introduction 
1.1 The Functional 

This paper is concerned with asymptotic properties of the following nonlocal Cahn-Hilliard 
energy functional defined on H 1 (R d ): 



1 

£ Jrpd 



£{u) := e J \Vu\ 2 dx + - / W(u)dx + 7 \\u — £ m||#-i(t<*)) (-'-•-'-) 



where we take the double- well potential W(u) := u 2 (l — u 2 ). Here the order parameter u 
is defined on the flat torus T d = R rf /Z rf , i.e. the square [—5, with periodic boundary 
conditions, and has two preferred states u = and u = 1. We are interested in the structure 
of minimizers of over u with fixed mass f Td u = f where / £ (0, 1). The first term e f |V«| 2 
penalizes large gradients, and acts as a counterbalance to the second term, smoothing the 
'interface' that separates the two phases. The third (nonlocal) term is defined as 



u 



IV-I'*, where - Aw — u - -f u. 

K ' JT d JT d 



This term favors high-frequency oscillation, as can be recognized in the l/|k| 2 -penalization 
in a Fourier representation: 

117/ lv\\ 2 ~ V |A(k)|2 

\\ u J u \\H- 1 (T d ) ~ 4vr 2| k |2- 

If the parameter 7 is large enough, this term may push the system away from large, bulky 
structures, and favor variation and oscillation at intermediate scales, i.e. give rise to patterns 
with an intrinsic length scale. As we explain in the sequel, we refer to this mass-constrained 
variational problem as the diblock copolymer problem. When the mass constraint / is close 
to or 1, minimizing patterns will consist of small inclusions of one phase in a large 'sea' of 
the other. We wish to explore this regime via the asymptotic behavior of the functional in a 
limit wherein 

• both e and the volume/mass fraction / of the minority phase tend to zero (appropriately 
slaved together) 

• 7 is chosen in order to keep the number of minority phase particles 0(1). 

We will primarily concern ourselves with the case d = 3 but remark on the analogous results 
for d = 2. 
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Figure 1: A two-dimensional cartoon of small-particle structures 



1.2 The Spherical Phase in Diblock Copolymers 

The functional £ was introduced by Ohta and Kawasaki to model self-assembly of diblock 
copolymers |23[ [22] . The nonlocal term is associated with long-range interactions and con- 
nectivity of the sub-chains in the diblock copolymer macromoleculd^ The order parameter 
u represents the relative monomer density, with u = corresponding to a pure-^4 region and 
u = 1 to a pure-i? region. The interpretation of / is therefore the relative abundance of the 
A-parts of the molecules, or equivalently the volume fraction of the ^-region. The constraint 
of fixed average / reflects that in an experiment the composition of the molecules is part of 
the preparation and does not change during the course of the experiment. From (1.1) the 



incentive for pattern formation is clear: the first term penalizes oscillation, the second term 
favors separation into regions of u = and u = 1, and the third favors rapid oscillation. Under 
the mass constraint the three can not vanish simultaneously, and the net effect is to set a fine 
scale structure depending on e, 7 and /. The precise geometry of the phase separation (i.e. 



the information contained in a minimizer of (1.1)) depends largely on the volume fraction /. 
In fact, as explained in [9], the two natural parameters controlling the phase diagram are 
r = (s 3 ^ 2 ^fy)~ 1 and /. When T is large and / is close to or 1, numerical experiments [9] 
and experimental observations [I] reveal structures resembling small well-separated spherical 
regions of the minority phase. We often refer to such small regions as particles, and they are 
the central objects of study of this paper. Since we are interested in a regime of small volume 
fraction, it seems natural to seek asymptotic results. Building on our previous work in [8], it 
is the purpose of this article to give a rigorous asymptotic description of the energy in a limit 
wherein the volume fraction tends to zero but where the number of particles in a minimizer 
remains 0(1). That is, we examine the limit where minimizers converge to weighted Dirac 
delta point measures and seek effective energetic descriptions for their positioning and local 
structure. 

The small particle structures of this paper are illustrated (for two space dimensions) in 
Figure [l] There are three length scales involved: the large scale of the periodic box T d , the 
intermediate scale of the droplets, and the smallest scale of the thickness of the interface. 
Two of these scales are known beforehand: the the size of the box we have chosen to be 1, 



1 See [10] for a derivation and the relationship to the physical material parameters and basic models for 
inhomogeneous polymers. Usually the wells are taken to be ±1 representing pure phases of A and B-rich 
regions. For convenience, we have rescaled to wells at and 1. 
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and the interfacial thickness should be 0(e) by the discussion above. The intermediate scale 
£, the size of the droplets, is not yet fixed, and will depend on the two remaining parameters: 
the parameter 7 in E and the volume fraction /. 

For a function u, the mass is defined as / = j Td u. In Fig. [I] the region where u ~ 1 
is small, suggesting that J Td u is small. We characterize this by introducing a parameter 77, 
which will tend to zero, and by assuming that the mass f Td u tends to zero at the rate of rj d : 



f 



M V d , 



for some fixed M > 0. 



(1.2) 



After rescaling with respect to rj, M will be the mass of the rescaled functions. We now have 
three parameters e, 7, and rj, which together determine the behavior of structures under the 
energy E etU ; in this paper we keep the parameter M fixed. Let us fix d = 3. In Section [3] we 
see that in terms of v := u/rf , the relevant functional is 



e rj 



\Vvfdx + 



T 3 



£ Jt 3 



W{v) dx 



+ r)\\v 



_ff _1 (T 3 ) ' 



where W(v ) := v 2 (l — rj 3 v) 2 . Via a suitable slaving of e to rj (see Theorem 3.1 ), we prove, via 
T-convergence, a rigorous asymptotic expansion for E^y^ of the form: 



E 



Eo + V^o + higher order terms, 



where both Eo and Fo are defined over weighted Dirac point masses and may be viewed as 
effective energies at the first and second order. Their essential properties can be summarized 
as follows: 

• Eo, the effective energy at the highest level, is entirely local: it is the sum of local 
energies of each particle, and is blind to the spatial distribution of the particles. The 
particle effective energy only depends on the mass of that particle. 

• Fo, the effective energy at the next level, contains a Coulomb-like interaction between 
the particles. It is this latter part of the energy which we expect enforces a periodic 
array of particles. 



The proof a Theorem [37L] relies heavily on our previous work for the sharp interface limiting 
functional E v (see Section [4] for its precise definition) obtained by fixing rj in E £:V and letting 
e tend to zero. The well-known Modica-Mortola Theorem [19] makes this limit E^ precise 
in the sense of T-convergence. The small-?/ asymptotics of E v were proved in [8j, and the 



main result of this article (Theorem 3.1) is to establish the same limiting behavior but in the 
diagonal limit of both e and r] tending to zero. We summarize these limits (for the leading 
order) in the diagram below. 

E, 

rj fixed, e— >0: 
Modica-Mortola Theorem 




Theorem|4.1| proved in [H] 



The article is organized as follows. In Section [3j we discuss the rescalings and state the main 
result Theorem 3.1 Section [4] explicitly states the main results of our previous paper [8] 



4 



which form the basis for the proof of Theorem 3.1 presented in Section [5} In Section [6j we 
discuss the variational problem associated with the first order T— limit Eo, connecting it with 
an old problem of Poincare and presenting some conjectures. In Section [7J we discuss the 
necessary modifications in two dimensions. 



2 Some definitions and notation 

We recall the definitions and notation of [8] . We use = H d /Zi d to denote the d-dimensional 
flat torus of unit volume. We will primarily be concerned with the case d = 3. For the use of 
convolution we note that T d is an additive group, with neutral element G T d (the 'origin' 
of T d ). For v G BV(T d ; {0, 1}) we denote by 

\Vv\ 

"pd 

the total variation measure evaluated on T d , i.e. ||Vu||(T d ) (see e.g. [2] or [31 Ch. 3]). Since 
v is the characteristic function of some set A, it is simply the notion of its perimeter. Let X 
denote the space of Radon measures on T d . For G X, ^ — v yu denotes weak-* measure 
convergence, i.e. 

/ fdfi v -)• / fdfi 

for all / G C(T n ). We use the same notation for functions, i.e. when writing — Vq, we 
interpret v„ and Vq as measures whenever necessary. 

We introduce the Green's function G T d for —A in dimension d on T d . It is the solution 

of 

— AG T d = 8-1, with / G T d = 0, 
where 5 is the Dirac delta function at the origin. In three dimensions^ we have 

G Ti (x) = + 5 (3) (x) (2.2) 

47r|x| 

for all x = (xi,X2,Xs) G R 3 with max{|xi|, \x2\, \%3\} < 1/2, where the function is 
continuous on [-1/2, 1/2] 3 and smooth in a neighborhood of the origin. 
For (i£l such that /^(T d ) = 0, we may solve 

— Au> = /JL, 

in the sense of distributions on T d . If w G H^T"*), then /i G tf-^T*), and 



ImIIh-ict"*) : = / \^w\ 2 dx. 



2 In two dimensions, the Green's function G T 2 satisfies 



for all x = (atijOia) £ R 2 with max{|a;i|, < 1/2, where the function t/' 2 ' is continuous on [—1/2, 1/2] 2 and 
C°° in a neighborhood of the origin. 
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In particular, if u G L 2 (T d ) then (u-fu) € H~ x (T d ) and 



L( T d) 



T »/ T 



u{x)u{y) G T d {x — y) dx dy. 



Note that on the right-hand side we may write the function u rather than its zero-average 
version u — -f u, since the function G T d itself is chosen to have zero average. 

If / is the characteristic function of a set of finite perimeter on all of R 3 , we define 



mm 



dx dy. 



/R3 JR3 47TF - V\ 

3 Rescalings and Statements of the Results 



We now rescale the energy £ in (1.1). Starting in three dimensions, for rj > we define 



v :- 



u 

^3 • 



so that £ becomes in terms of v 
e rj 

where 



Vv\ 2 dx + — 



T 3 



e JT 3 



W(v)dx + 7?? 6 ||v - f u[|fl--i( T 3), 



(3.1) 



W{v) := v 2 (l-rfv) 2 



In order to find the correct scaling of 7 in terms of 77, we consider a collection u„ : T 3 



{0, 1/rj } of components of the form 



(3.2) 



where the are disjoint, connected subsets of T 3 . Then under the assumption that the 
number of A4 remains 0(1), we find 



e 77 



Nvrdx + 



T. 



W(v) dx 



~ 77 



0(1). 



Here we are using the well-known Modica-Mortola convergence theorem \19\ [5] linking the 
perimeter to the scaled Cahn-Hilliard terms. A simple calculation (done in [8]) shows that 



the leading order of the ||v„ — f 1 



'r)llff-l(T3) 



is I/77, and that this leading contribution is from 



the self -interactions, i.e. — j- ^||#-i( T 3) is I/ 7 /- Thus balancing the third term in (3.1) 
implies choosing 7 ~ I/77 3 . Hence we set 



77 



3 ' 



Choosing the proportionality constant equal to 1 entails no loss of generality, since in the 
limit e — > this constant can be scaled into the mass M defined in (1.2). 
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With this choice, one finds 



£(u) = rj 1 < 77 



\Vv\ 2 dx + — 



T 3 



e JT 3 



W(v) dx 



(T3) 



noting that the contents of the outer parentheses is 0(1) as r\ — > with e <^ r/. Thus let us 
define the re-normalized energy 

,3 



E e , v (v) := rj 



£ rj 



\Vv\ 2 dx + — 



T 3 



£ JT 3 



L (T 3). 



(3.3) 



We are interested in the small-77 behavior of E £ ^ and describe this behavior via functionals 
defined over Dirac point masses. Let us first introduce the remaining relevant functionals in 
our analysis. First we define the surface tension 



a := 2 / y/W(t) dt. 
Jo 



For the leading order, we define 

eo(ro) :=inf j<rjfjVz| + IM|^-i (R3) : z G W(R 3 ; {0, 1}), J 



z = m 



and 



Eo(v) := 



R 3 



e o(^ J ) if f = m% ^x i i { x% } distinct, and m l > 
00 otherwise. 



(3.4) 



(3.5) 



For the next order we note that among all measures of mass M, the global minimum of Eo is 
given by 



min I E Q (v) : j v = M j = e (M). 



We will recover the next term in the expansion as the limit of E ert — eo, appropriately rescaled, 
that is of the functional 



Its limiting behavior will be characterized by the functional 

- 00 

J> (3) (o) (m i ) 2 + 



Fo(«) := 4 



i=l 



m 1 ™? G T z(x l — xi) if f = m l 5 x i with {x 1 } distinct, {m 1 } £ 

i=l 

otherwise, 



where g( 3 ) is defined in (2.2) and 
M := 



{m ? }igN : m l > 0, eo^ 4 ) admits a minimizer for each i, 

00 

and ^^eo(w- i ) = e o(^2 



m 



i=l 



i=l 
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In Theorem 3.1 we prove that 



r r 
E £>r) — > E Q and F £>ri — > F . 



Precisely, 



Theorem 3.1. • (Condition 1 - the lower bound and compactness) Let e n and r/ n be 

sequences tending to zero such that, for some ( > 0, e n = o(r/t +< *). Let v n be a sequence 
such that the energy E En>r)n (v n ) is bounded. Then (up to a subsequence) v n — vq, suppt>o 
is countable, and 

liminf E £ Vn (v n ) > E (v ). (3.6) 

n— yoo 

If in addition E £n ^ n (v n ) is bounded and £ > 1, then the limit vq is a global minimizer 
of Eo under constrained mass, and 

liminf F £ntVn (v n ) > Fq(v ). (3.7) 



• (Condition 2 - the upper bound) There exist two continuous functions C\,Ci : [0, oo) — > 
[0, oo) with Ci(0) = 6*2(0) = with the following property. Let e n and rj n be sequences 
tending to zero and let e n < Ci(r] n ). Let vq be such that Eq(vq) < 00. Then there exists 
a sequence v n — vq such that 

limsup^^^n) < E («o)- (3.8) 

n— >oo 

If in addition vq minimizes Eq under constrained mass and e n < C^r/n); then this 
sequence also satisfies 

lirasupF £ntVn (v n ) <F (v ). (3.9) 

n— >oo 



Remark 3.2. Choice of the slaving of e to r/. There are two separate arguments 
connecting the two parameters: 

• If the sharp-interface approximation is to be reasonable, then the scaling should be 
such that the interfacial width is small with respect to the size of the particles. Since a 
particle has diameter 0(r/), this translates into the condition e -C rj. 

• Eo is infinite on structures that are not collections of point masses. If Eo is to be the limit 
functional of E E „, then along any sequence that does not converge to such point-mass 
structures E £)Tj should diverge. It turns out that this provides a stronger condition, as 
we now show. 

For E er) , every function v is admissible. Under constrained mass M, an obvious candi- 
date for the limit behavior is the function v = M, with energy scaling E er ,(X) ~ ?? 4 /e. On 
the other hand, if the functional E £tV is close to Eo, then we will have E £<t] ~ Eo = 0(1)- 
Therefore the ratio rj^/e is critical. If this ratio is small, then the constant state has 
lower energy than localized states, and we do not expect the functional E v to be a good 
approximation of E £r) . On the other hand, if the ratio rj A /e is large, then localized 
states have lower energy than constant states. 
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In Theorem 3.1 above, the lower bound is responsible for forcing divergence of the energy 
along sequences which do not converge to point masses; the lower bound therefore requires 
£ "C rf . The extra factor r/£ is used in the truncation part of proof: in relating a diffuse- 
interface sequence to a sharp-interface sequence, we truncate at a suitable level set of the 
interface, and the small factor rjn is used to quantify the closeness in interfacial energies with 
respect to the surface tension a. 

For the upper bound, we would ideally require e n = o(%). What we assume, e n < Ci(j] n ) 
and e n < 62(77™), are stronger requirements. However, at this stage we do not know the exact 
local behavior for minimizers of eo- In two dimensions we can fully characterize this local 
behavior, and as we shall see in Section [7j this allows us to require only the weaker condition 
£n = o(rj n ) (up to a logarithmic correction). In three dimensions we use a convenient version 
of the Modica-Mortola profile construction which does not give an optimal scaling in terms 
of closeness of energies (cf. Lemma 5.1 ). If one can establish the conjectured behavior for the 
local problem (see Section [6]), one can then achieve a sharper slaving. 



4 Previous results for the sharp interface limit 

In [8] we dealt with the sharp-interface functionals that arise from letting e tend to zero for 
fixed n. For E £ ^ and F E „ respectively these limit functionals are 



E, := 



V/t3I V "I + v\\v-fvf H _ 1{T3) if v e BV(T 3 ; {0, 1/rr 3 }) 
00 otherwise, 



and 



E v (v) - e v) 



We proved that 



r r 
E v — > Eo and — > Fq, as 77 — >• 0. 



Precisely, 

Theorem 4.1. Let r\ n be a sequence tending to 0. 

• (Condition 1 - the lower bound and compactness) Let v n be a sequence such that the 
sequence of energies E Vn (v n ) is bounded. Then (up to a subsequence) v n — ^ vq, supp^o 
is countable, and 

liminf E Vn (v n ) > Eo(«o). (4.2) 

n— >oo 

If in addition E Vn (v n ) is bounded, then the limit vq is a global minimizer of Eo under 
constrained mass, vq = ^2ini l 5 Xi where {m ! } £ A4 and 

liminf F r?n (?j n ) > F (f )- (4.3) 

n— >-oo 

• (Condition 2 - the upper bound) Let Eo(^o) < 00 and Eq(vq) < 00 respectively. Then 
there exists a sequence v n — vq such that 

limsupE^ n (7; n ) < E («o)- (4.4) 
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If Fo(wo) < 00 then exists a sequence v n vq such that 

limsupF^ n (t; n ) < F (vo)- (4.5) 

We recall from [8] some properties of e§: 

1. For every a > 0, e is non- negative and bounded from above on [a, 00). 

2. If {m*}j g N with Ylii mt < 00 satisfies 

00 00 

^e (m J ) = eo (^m i ), (4.6) 

i=l i=l 

then only a finite number of m l are non-zero. 



Remark 4.2. In proving Theorem 4.1, the bulk of the work was confined to the lower- 
bound inequalities wherein, after establishing compactness, one needed a characterization of 
sequences with bounded energy and mass. The characterization implied that such a sequence 
eventually consists of collection of non-overlapping, well-separated connected components 
(see [U Lemma 5.2]). 

We note that in proving the second-order V convergence, we saw that for an admissible 
sequence v n , the boundedness of F Vn (v n ) implied both a minimality condition and compact- 
ness: 

• The minimality condition arose from the fact that E £n ^ n (v n ) must converge to its 



minimal value and implied that the {m ! } must satisfy (4.6). Hence by property 2 
above, the number of limiting particles must be finite. 

The compactness condition implied that for each m l , the minimization problem defining 



eQ(m % ) (namely (3.5)) had a solution. 



These condition are responsible for the additional properties of the weights m' (c.f. Ai) in 
the definition of Fq. 



5 Proof of Theorem 3.1 



The proof of Theorem 3.1 relies on Theorem 4.1 For the lower bound we use a suitable 
truncation to relate the approximating diffuse-interface sequence to a sharp-interface sequence 
with the same limit and whose difference in energy is small. For the upper bound, we modify, 
in a neighborhood of the boundary, the sharp-interface recovery sequence given by Theorem 

Such a 



4.1 via a quantification of the Modica-Mortola optimal-profile construction ( [19 
result is provided by a lemma of Otto and Viehmann [24J. 

Lemma 5.1. Let a > 0. There exists a constant Co (a) such that for any characteristic 
function x of a subset of T 3 and 5 > 0, there exists an approximation u G ff 1 (T", [0, 1]) with 

1 



8\Vu\' + 



■ u 



(1 



and 



I 



\x — u\ dx < 



1 dx < (a + a) 



C Q (a)5 I |V X |. 
Jt 3 



T :1 



|Vx|, 
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The proof of Lemma 5.1 follows from the proof of Proposition 1, Section 7 in [25]. Note 
that in [23], the authors deal with the functional 



n2(l-u 2 ) 



IVttl 



+ 2S (1 



u 2 ) dx, 



defined on cubes of arbitrary size Here the wells are at ±1 and more importantly, this 
scaling produces unity as the limiting surface tension a. However the structure of their proof 
only uses the fact that this functional T-converges to 



iVd. 



Hence our Lemma 5.1 follows directly not from the statement of their Proposition 1 but from 
its proof. 



Proof of Theorem 3.1, We first prove Condition 1 (the compactness and lower bounds). Let 
e n , rj n , and v n be sequences as in the theorem, such that E £n ^ n (v ri ) is bounded (but not 
necessarily F £n>rin (v n ), yet). For part of the proof we will work with the sequence and the 
energy in the original scaling u n , given by u n = Vn v n- In terms of u n , we find 

E £ n ,Vn( V n) = % I |Vn n | 2 + ^ 1 — / W(u n ) + -n\\u n - j Mnllff-l, 
Vn Jt* Vfen J T 3 Vn 

Following ; 1T)J we define the continuous and strictly increasing function 

and note that as a consequence of the inequality a 2 + b 2 > 2ab, we have 

If 1 

E £n , Vn (v n ) > ^ \V(f)(u n )\ + -^\\u n - funW 2 ^. (5.1) 

Vn JT 3 Vn 

Now set a n = l/(a — Vn), where as before a = 2 ^/W(t) dt = 4>(1) — 0(0). Fix 5 n > 
by the condition 

0(1 - 26 n ) - 0(2<5 n ) = 0(1) - 0(0) -vi = ~, 

a n 

and note that the quadratic behavior of W at and 1 implies that 5 n = 0(rj n )• We also 
introduce the notation [it] for the clipping to the interval [0, 1]: 

[it] := min{l, max{0, it}}. 

Using the characterization of perimeter (cf. [12] or [2j Th. 2.1]) as 

/ |V0(W)| = ^ U\d*{4>{[Un]) > t})dt, 

Jt 3 J 4,(0) 

we estimate the size of the set 

A n := |t€ [0(0), 0(1)] :'H 1 (d*{HK])>t}) > a n ^ | V0([i/ n ] 
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by 

14.1=/ ldt< I f 1 H 1 (d*{<f>([u n ])>t})dt=—. 

JA n «n J T 3 |V0(M)| J 4,(0) OL n 

By the definition of a n and S n it follows that there exists a t n £ [4>(8 n ),4>(l — S n )] \ A n , for 
which therefore 

U\d*{ct>{[u n ]) > t n }) <a n f \V<f>([u n })\ <a n f \V<f>(u n )\. (5.2) 

We now construct an auxiliary sequence u n such that the corresponding v n = u n /r] n will 
be admissible for the sharp- interface functional E v . We map the values of u n to {0, 1} with 
cutoff 4>~ l (t n ): 

if (j)(u n (x)) < t n 

1 if <p(u n (x)) > t n , 



u n (x) : 

so that 



Vu n | =n 1 (d*{<P([u n ]) >t n }). (5.3) 



We estimate the difference in L 2 and H^ 1 of u n and u n . Since ^ -1 (i„) G [<5 n , 1 — <5 n ], the 
function 

1i 2 if 0(ti) < t n 

a-u) 2 if <f>{u)>t n 

is bounded from above by an increasing factor times W, 



M u ) - 



M u ) < C5- 2 W{u) < C'r]n C W{u) for some C, C' independent of n. 
Therefore the sequences u n and u n are close in L 2 : 

[ MO < C 'Vn C [ W{u n ) = 0(e n r] 2 n - C ) -> 0, 

JT 3 JT'-i 



I — II 2 

\u n — u n 1 1 ip. 



where the final estimate results from the boundedness of E £njVn (v n ). Consequently they are 
also close in H -1 , 

\\Un-U n -f(Un-Un)\\ H -i < C\\u n ~ U n ~ f{u n - U n ) \\ L 2 

< C\\u n -Un\\ L 2 

= 0{eHW-^) 0, (5.4) 

and the same holds for the squared norms: 

In r |i2 _ II— C — II 2 | 

\\\Un J^nufj- 1 W^n J^"n\\}{-i\ 

< (\\ u n -fu n \\ H -l + \\u n ~ f U n \\ H -l) \\U n -U n ~ j{u n - « n )||if-l 

< (2||n n -ju n \\ H -i + \\u n -u n - j{u n - « n )||if-i) 0{e 1 J 2 r] n ~ <i/2 ) 
= (v 5 n E £niT]n (v n )) 1/2 0(el/ 2 Vn -« 2 ) + 0(s nV 2 n -<) 

= 0(el/W n /2 - C/2 ) + 0(e nV l-<) 

= o(r, n ), (5.5) 
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Note that in the last lines of (5.4) and (5.5), we have used the hypothesis e n = o(rjn + ^). 



Using (5.2) and ( |5.3[ ) we transfer the lower bound (5.1) to the sequence u n : 



E, 



£n,rjn \ v n 



1 n\d*{^{[u n \) > t n }) + 4lK -fu n f H -i 

Vn 



( [5T3]t , ([5T5] ) 



a nVn 
1 



Vr. 
a, 



i« 1 1 „ 

V«n + -F \\Ur, 

Vn 



ju n \\ 2 H -i + o(r]n 



/ \Vv n \ + Vn\\ v n ~ j~ VnW'fj-i + °{Vn) 



> 



era, 



- E r?„(Wn) + o(T] n ), 



(5.6) 



where in the last line we used the fact that era n > 1 (note that era n —> 1 as n — > 00) 



From (5.6) it follows that the sequence v n satisfies the conditions of Theorem 4.1 There- 
fore there exists a subsequence v nk converging to a limit vo, with countable support, such 
that 



liminf E v (v nk ) > E (v ). 



k— >oc 



(5.7) 



The corresponding subsequence v nk of the sequence v n also converges weakly to the same 
limit, since for tp £ C(T 3 ), 



T 3 



1 



lfe || i2 ||^|| i2 = 0(4/2^/2) 



This proves the compactness of the sequence v n and the characterization of the support of 
the limit vq. The lower-bound inequality (3.6) then follows from (5.6) and ( |5.7| ). 

We address the lower bound for F £>rj . We note that boundedness of F £n ^ n (v n ) implies 
boundedness of E £n!Vn (v n ), so that the characterization of the convergence of the sequence 
given above applies. In addition, by (5.6), we have 



> 



1 

Vn 
1 

Vn 



Ee n , Vn (v n ) - e 



X 3 



+ 



1 



1 



Vn \cra. 



1 E,„(S„) + o(l). 



Since era n = 1 + o(r] n ), with ( > 1, the lower bound (4.3) for F„ implies 



lim inf F, 



(v n ) > Fo(^o), 



which is (3.7). 



We now turn to the upper bound (condition 2), treating E £t „ first. As in the proof of 
Theorem 4.1 it is sufficient to prove that for any vq of the form 

N 

vq = ^m'^i, with x 1 distinct, 
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there exists a sequence v n vq with 

lim sup J3 e „ < Eo(v ). 



(5.8) 



See [S] for an explanation. Given such a vq, Theorem |4.1| (specifically (4.4)) provides an 
admissible sequence v n — 1 Vq for with 



lim E Vn (v n ) = E (v ). 



(5.9) 



We write u n := Vn v ni which is the characteristic function of a subset of T 3 composed of N 
sets whose diameters are decreasing to zero. For each n, Lemma 5.1 with a = r/ n implies that 
there exists a Co(r) n ) such that for any e n > 0, we have an approximation u n G H X (T 3 , [0, 1]) 
such that 

/ e n |Vu n | 2 + — u 2 n {\ -u 2 n )dx < (a + r/n) / \Vu n \, 

JT 3 £ n Jt 3 



(5.10) 



and 



Now let 



We have 



T 3 



I dx < C (rj n )e n 



iVitJ. 



T 3 



\V n — UnUL 1 ^ 3 ) 



< 



Un 
Vn 



1 

Vn Jt 3 
Co(r?„)e„ 



Vn 



'Ufi I dec 
|V« n | 

T 3 



< c 



C (rin)e n 

Vn 



(5.11) 



We will slave e n to r\ n such that the above tends to zero as n tends to infinity. In particular, 
v n and v n will have the same limit Vq. We crudely estimate the i7 _1 -norm as follows 



\v n - v. 



n f( v n v ™)||_H"-i(T 3 ) 



< C\\Vr, 



W( 



V n - V n 



lL 2 (T 3 ) 



< C \\V n — V n \\ L 2(r£3j 

< C\\v n — U n ||i°°(T 3 ) \\ v n -V n \\ L i( T :^ 
^ c Coi^En 



Vn 



(5.12) 
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Next we note that 

E e n , Vn {Vn) 



< 



[ \Vu n \ 2 + ^— [ W(u n ) + \\\u n - j u n \\ 2 H -i 
T / (e n \Vu n \ 2 + — u 2 {\ - u 2 n ) ) dx + r] n \\v n - § v n \\ 2 H - 
~2 / [e n \Vu n \ 2 + — u 2 (l -u 2 )) dx + rj n \\v n - j v n \\ 2 H - 



v 
j 

Vn 
1 

v, 



([5TT0JJ5TT2J 



Co(?7n)£n 
»7n 



?7n(0" + ^n) / |V« n | + %||u n - f UnHjy-l + C 
JT3 



Thus we assume 



Er,M + ?7n / |V«„| + C 
Jt 3 

Co(^n)£n 



vl 



— >• as n — > 00, 



(5.13) 



(5.14) 



and we choose a function C\ as in the Theorem such that (5.14) is satisfied whenever e n < 
Ci(rj n ). We now take the limsup as n — ^ 00 in (5.13), and hence (5.9) gives (5.8). 

For the next order, let 

v 



i=l 



Theorem 4.1 (specifically (4.5)) gives a sequence v n — i>o such that 

lim F Vn (v n ) = F {v ). 



(5.15) 



We take v n to be the diffuse-interface approximation used in the previous upper-bound ar- 
gument but now taking a to be rj 2 . Hence v n — 1 vo and following the steps of (5.13), we 
have 

Co(Vn) £ n 



E en , Vn {vn) < E„K) + rf n [ \Vv n \ + C 

JT 3 Vn 



(5.16) 



and 



E„K) + ^ / |V^| + C 

JT3 Vn 



eo 



T 3 



v n + e 



T 3 



eo 



T 3 



< F^K) + 0(r/ n ) + r/" 1 



£ \\v n ~ V n \\ L l + C 



Co(vl)er 

Vn 



where L is the local Lipschitz constant of eo (cf. Remark 4.2). Thus choosing e n such that 

Co(vl)£n 



Vn 



as n — > 00, 



(5.17) 
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(5.15) implies 



limsupF £n >Vn {v n ) < Fq(v ). 



We choose a function C2 as in the Theorem such that e n < C2{rj n ) implies (5.17). 



□ 



6 The local structure of minimizers and the variational prob- 
lem that defines eo 

Simulations of minimizers of the diblock copolymer problem show phase boundaries which 
resemble constant mean curvature surfaces (see for example [9] and the references therein): 
In the regime of this article, we observe spherical boundaries. Experimental observations in 
diblock copolymer melts also support this |32j . On the other hand one can see, for example 
via vanishing first variation, that on a finite domain the nonlocal term will have an effect 
on the structure of the phase boundary |llj . While rigorous results on this effect remain 
open, one would expect that exploiting a small parameter might prove useful and, indeed, 
this is exactly what our first order asymptotics have done: in proving the first order lower 
bound, we have reduced the local optimal shape of the particles to solutions of the variational 



problem (3.5) that defines eo- The details of this calculations can be found in (8]. Let us now 



comment of this problem and present some conjectures. 

We briefly recall the problem defining eo- For m > 0, minimize 



\Vu\ + 



u(x)u(y) 
4ir\x — y\ 



dx dy over all u e BV(R 3 , {0,1}) with / udx 



rn. 



Note that the two terms are in direct competition: balls are best for the first term and worst 
for the seconcj^J The function eo(m) denotes this minimal value, i.e. 



eo(ra) := inf < / |Vu| + 



u(x)u{y) 
4tt\x — y\ 



dx dy 



u G W(M 3 ,{0,1}), / udx 



m 



The latter point has an interesting history. Poincare [251 126] considered the problem of determining 
possible shapes of a fluid body of mass m in equilibrium. Assuming vanishing total angular momentum, the 
total potential energy in terms of it, the characteristic function of the body, is given by 



CP) 



u{x) u(y) 



C\x- 



dx dy, 



where — (C\x — y\) , C > is the potential resulting from the gravitational attraction between two points x 
and y in the fluid. Poincare showed under some smoothness assumptions that a body has the lowest energy if 
and only if it is a ball. He referred to some previous work of Liapunoff but was critical of its incompleteness. 
It was not until almost a century later that the essential details were sorted out wherein the heart of proving 
the statement lies in rearrangement ideas of Steiner for the isoperimetric inequality. These ideas are captured 
in the Riesz Rearrangement Inequality and its development (c.f. [H]): for functions /, g and h defined on R d , 



f{y)g{x-y)h{x)dydx < 



fr 



{y)g*{x~y) h*(x)dydx, 



where /*,<?*, h* denote the spherically decreasing rearrangements of /, g, h. While the general case of equality 
was treated by Burchard in [7], for the problem at hand where the function g ~ ■ j _1 is fixed and symmetrically 
decreasing, the inequality with the specific case of equality was treated by Lieb in [T7], thus proving that balls 
are the unique minimizers for the potential problem (V). 
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We also define the energy of one ball of volume m: 



2/3 



2/3 



/(to) := (Seyr) 1 / 3 ™ 2 / 3 + - to 5 / 3 . 

5 \47ry 

Clearly, we have eo(m) < f(m). We conjecture the following scenario. There exists to* > 0, 
such that for all m < to*, there exists a global minimizer associated with eo(m), and it is a 
single ball of mass m. For m > m*, a minimizer fails to exist. In fact, as to increases past 
to*, the ball remains a local minimizer, but a minimizing sequence consisting of two balls 
of equal size that move away from each other has lower limiting energy. This separation is 
driven by the H -1 interaction energy, which attaches a positive penalty to any two objects 
at finite distance from each other. The limiting energy of such a sequence is simply the sum 
of the energies of two non-interacting balls, i.e. 2/(m/2). The critical m* then is the only 
positive zero of f(m) — 2/(m/2), m* ~ 22.066. 

As m further increases above a certain m** > to*, a sequence consisting of three balls 
of equal size is a minimizing sequence for eo(m), with limiting value 3/(m/3); and so on for 
higher values of n. Specifically, we conjecture the following: 

Conjecture 6.1. The minimizer associated with eo(m) exists iff to < to*, and it is a ball of 
mass to. Moreover, for all m > 0, we have 

eo(m) = inf nf(m/n). 

ngN 

The infimum is achieved iff to < to*. 

One might ask as to what is known about global minimizers. In our previous article [8] on 
the sharp-interface functionals, we prove that if a sequence (in n) has bounded energy F v then 
it must converge to a weighted sum of delta functions where all the weights m l must have 
a corresponding minimizer of eo (to*) . One can readily check, via trial functions, that such a 
sequence exists. Thus for certain values of to, a minimizer of eo(rn) does exist. Unfortunately, 
our lower bound compactness argument gives no explicit range for the possible limiting weights 
to*. One could also consider local minimizers, and in particular one can study the stability of 
balls. A calculation of the second variation using spherical harmonics (unpublished) indicates 
that the ball retains stability up to m c ~ 62.83, well past the critical mass to*. 



Proving Conjecture 6.1 would for the first time provide some rigorous justification for 
why minimizers of the diblock copolymer problem have phase boundaries which resemble 
periodic constant mean curvature surfaces, supporting the idea that at small length scales, 
the perimeter (short-range) effects override the nonlocal (long-range) effects. 



7 Analogous results in two dimensions 

As in [8j, we summarize the analogous results for d = 2. While we do not give all the details, 
we give the essential features which should enable the reader to complete the proofs. The 
fundamental difference between two and three dimensions is that the if _1 -norm is critical in 
two dimensions. As explained in [8j, after rescaling with v = u/rj 2 this involves slaving 7 to r\ 
via 

1 

^ I log 771 77 s ' 
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I .fr- 



aud the two-dimensional function analogous to E E „ becomes 

E™(v) :=erf j \Vv\ 2 + ^ j W(v) + |log r/p 1 \\v - f • 

Here the rescaled double-well energy is now 

W{v) := v 2 (l-r] 2 v) 2 . 
The analogous sharp-interface (e — > 0) limit is given by 

ariJ T \Vv\ + llogr?!" 1 \\v - j v\\ 2 H _ 1(T) if v € W(T, {0, 1/r? 2 }) 



oo 



otherwise, 



where a is again given by (3.4). The first-order limit is defined by 

Yliel eo D ( m *) if ^ = Yl l i&i m% ^x i with I countable, {x % } distinct, and m i > 



Ef (v) :- 



x 



otherwise. 



where the function : [0, oo) — > [0, oo) is defined as follows. Let 

,2 



m 



e n (m) := h 2a\/irm 

47T 



777- 

h inf <! a I \Vz\ : z £ 

4^ 1 7 R 2 



W(R 2 ;{0,1}), / 



z = m 



(7.1) 



An interesting feature here is the explicit nature of e 2 ^ (in contrast to (3.5)). The first 
term is the dominant part of the H^ 1 norm in two dimensions, and it arises from the fact 
that the logarithm is additive with respect to multiplicative scaling. We introduce the lower- 
semicontinuous envelope function 



eg°(m) : = inf <^ eg D (m J ) : m? > 0, 



m J = m 



(7.2) 



j6J 



For the next order, note that 



min<! ElPiv) : 



v = M 



T 2 



e™(M). 



We hence recover the next term in the expansion as the limit of E 2D — efP, appropriately 
rescaled, that is of the functional 



Note that the corresponding sharp interface function is 



Ff 0) : = llogr/l 



In order to define the second-order limit, we require some preliminary definitions. We first 
recall a lemma whose proof was presented in [8] . 
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Lemma 7.1. Let {m 2 }jgN be a solution of the minimization problem 



mm 



j^eg D (m < ) : > 0, 



. i=l 



oo 

E 

i=i 



m 



M. 



(7.3) 



Then only a finite number of the terms m l are non-zero and all the non-zero terms are equal. 
In addition, if one m l is less than 2~ 2 / 3 7r, then it is the only non-zero term. 



Let 

For n£N and m > the sequence n < 

(n <g) m) 4 



m 



3 - 2 log 
m is defined by 



m 



1 < i < n 

n + 1 < % < oo. 



Let Ai be the set of optimal sequences for the problem (7.3): 

M. : = |n ®m:n®m minimizes (7.3) for M = nm, and eQ D (m) = e 
Then define 

'n{/ (m) + mV 2 )(0)} + 



2D, 



m 

T 



X J 



i,j>l 
i+3 



if v = m &x i j {x 1 } distinct, n®m G jM, (7 4) 



1=1 

otherwise, 



where the function g( 2 ) was defined in (2.1). We briefly comment on these functionals and 
their properties. As in three dimensions, the boundedness of implies that the limiting 
weights m 1 satisfy both a minimality condition and a compactness condition. The minimality 
condition implies that n ® m minimizes (|7.3l) . The compactness condition implies that 



,2D, 



m 1 ) 



e 2D K 



(7.5) 



As we can see from Lemma 7. 1 , the minimality condition provides a characterization that 
is stronger than in three dimensions: in particular the masses must be equal. Let us also 
comment on the function /o- The minimization problem (7.1) has only balls (here circular 
disks) as solutions. Thus in computing the small-?7 asymptotics of F^, the i/ _1 (R 2 )-norm 
of a two-dimensional disc of mass m enters. The functional fo(m) is exactly this value. 



Theorem 7.2. 



(Condition 1 - the lower bound and compactness) Let e n and rj n be 



sequences tending to zero such that £ n r/, 



-3-C 



for some £ > 0. Let v n be a sequence 



such that the energy E^„ n (v n ) is bounded. Then (up to a subsequence) v n — vq, supp vo 
is countable, and 

liminf E™ {v n ) > Ef{v G ). (7.6) 

n^oo ' ' 

// in addition F^„ (v n ) is bounded, then the limit Vq is a global minimizer of Eg D under 
constrained mass, and 

liminf F^M > F§ D («„). (7.7) 
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( Condition 2 - the upper bound) Let e n and r\ n be sequences tending to zero such that 
e-nVn 1 \1°&Vn\ ~ > 0- Let v be such that Eq B (v) < oo. Then there exists a sequence v n — v 
such that 

n— >oo 

If in addition v minimizes Eg D under constrained mass, and if Erj}~ x |log?7 n | 2 — > 0, then 
this sequence also satisfies 



The proof of Theorem 7.2 is very similar to that of Theorem 3.1 Again, we rely heavily on 
the lower bound estimate and upper bound recovery sequence of the associate sharp interface 
problems. We summarized those results in [8]. The lower bound inequality follows verbatim 
the three-dimensional case, the differences in dimension reflected by the exponent 3 as opposed 
to 4 in the slaving of e n to rj n . 

The main difference comes in the upper bound and this is reflected in the less restrictive 
slaving of e n to rj n . In two dimensions, minimizers associated with the first order limit are 
necessarily circular droplets. This gives an upper-bound recovery sequence of circular droplets 



(cf. (7.1)). To regularize the circular boundaries, one can bypass Lemma 5.1 and simply use 



a one-dimensional optimal profile to approximate a Heaviside function. The advantage here 



is the explicit dependence of Co on r\ n . For the analogous step to (5.12), one can use the 
following interpolation inequality corresponding to the 'nearly'-embedding of L 1 in H^ 1 to 
relate the i? _1 -norm to the L 1 -norm: 

Lemma 7.3. Let f 6 L°°(T 2 ) with f T2 f = 0. Then there exists a constant C > such that 

\\f\\l- w <c\\f\\h iT2) (i + ^ ||J||L1(T2)/ 

Since the proof this inequality is short and to our knowledge, absent from the literature, 
we end this section with its proof. To this end, we first derive an inequality proved by Brezis 
and Merle [6] in a slightly different form. 



Lemma 7.4. There exists a constant Cq > 1 such that 

[ e^<C , 

for all (j) £ VF 2,1 (T 2 ) satisfying 



|A</)| = 1. 

T 2 

Remark 7.5. As the proof below shows, the result holds true for any (j> such that J" T2 | A0| < 
47r; the constant Cq diverges as the critical value of 4-zr is approached. 

Proof of Lemma 7.4 Setting f(x) := — A(f>, so that J \ f\ = 1, we have 

(p(x)= G T 2(y)f(x - y)dy, 
Jt 2 
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and note that by (2.1) 

|G T3 (y)| < C- ^log|y| 
for some C > 0, and for all y £ (-1/2, 1/2) 2 . Therefore, using Jensen's inequality, 

/ e^dx< I expf/ \G T2 (y)\\f(x-y)\dy ) dx 

JT 2 J(-l/2, 1/2) 2 W(-l/2,l/2) 2 / 



log(|yr 1/2 ") \f(x-y)\dy) dx 

-1/2,1/2)2 V / / 



<e c exp 

J (-1/2,1/2)* \J (-1/2,1/2) 

<e c [ [ \y\-V 2 *\f(x-y)\dydx 

J (-1/2, 1/2)2 J(_i/ 2j i/2)2 

= e c [ \y\- 1/2n dx 

J(-l/2,l/2)2 

= :C . 



□ 



Proof of Lemma \7.3[ Set 3>(s) := |s| log(l+Cb|s|) and let 3>* be the convex conjugate $*(i) := 
sup sgR (ts — 3>(s)). From the lower bound <3?(s) > |s| log(Cb|s|) we derive the upper bound 

$*(t) < Cq-^I*!. 

Define the Orlicz norm 

inf <! A > : / (-) < 1 



_> 



.A/ 

Then we have the Holder inequality (see, for example, Section 3.3 of [27J 



/ fg<2\\fh\\gU*- 

Jt 2 



To prove Lemma 7.3 we take / £ L°°, / ^ 0, with J / = 0, and by multiplying / with a 
constant we can assume that J \ f\ = 1. Setting —Ac/) = /, we have 



2 

H-!(T2) 



/0<2||/||*|M|*.< 2||/||*. 

T 2 



The second inequality above follows from remarking that 

>. =inf{A>0: [ $*(|) <1 

<mf|A>0:C ( T i y e^ l/A < 1 

LemmaJ73] 
< 1. 

Now let A* := ||/||$. Since the map A — > J T - 2 ^(a) * s continuous at A*, we must have 

X 2 V A* 
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Thus 



\* = \f(x)\ log + C ^ j d.v £ log ( 1 + ( n 



or 



fix 

A* 

A*(e A * - 1) < C 



We note that 

1<>S2 e A < A(e A - 1) for all A > with A(e A - 1) > 1. 



2 

Hence if A*(e A * - 1) > 1, then 

||/|| 2 H- 1(T2) < 2||/||. = 2A, < 21og 2 ^ g2 , 

On the other hand, if A*(e A * — 1) < 1, then since A i— > A (e A — 1) is increasing, we have A* < A, 
where A(e A — 1) = 1. Since 

e*<A(e*-l), 

we have 



&_ 1(T3) < 2||/||* = 2A* < 2A < 2 log 
Replacing / with //||/||x,i gives the desired inequality. 



□ 



8 Discussion, dynamics, and related work 

Together with [Sj, we have presented an analysis of the small- volume regime for the diblock 
copolymer problem. This has been accomplished by an asymptotic description of the energy 
functional in the small volume- fraction regime. We refer to the discussion section of [8] for 
comments on the role of the mass constraint with respect to the limit functionals and the 
fundamental differences between the two- and three-dimensional cases. As described above, 
in three dimensions many open problems remain with respect to the local structure problem 
and it is here that one should first focus in order to rigorously address the role of the nonlocal 
term on shape effects. 

This asymptotic study has much in common with the asymptotic analysis of the well- 
known Ginzburg-Landau functional for the study of magnetic vortices (c/. |3CH [15l Q]). Our 
problem is much more direct as it pertains to the asymptotics of the support of minimiz- 
ers. This is in strong contrast to the Ginzburg-Landau functional wherein one is concerned 
with an intrinsic vorticity quantity which is captured via a certain gauge-invariant Jacobian 
determinant of the order parameter. 

Our results are consistent with and complementary to two other recent studies in the 
regime of small volume fraction. In [28] Ren and Wei prove the existence of sphere-like 



solutions to the Euler-Lagrange equation of ( 1.1 ), and further investigate their stability. They 



also show that the centers of sphere-like solutions are close to global minimizers of an effective 
energy defined over delta measures which includes both a local energy defined over each point 
measure, and a Green's function interaction term which sets their location. While their results 
are similar in spirit to ours, they are based upon completely different techniques which are 
local rather than global. Very recently. Muratov [20] proved a strong and rather striking 
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result for the sharp interface problem in two dimensions. In an analogous small volume 
fraction regime, he proves that the global minimizers are nearly identical circular droplets of 
a small size separated by large distances. While this result does not precisely determine the 
placement of the droplets - ideally proving periodicity of the ground state, to our knowledge it 
presents the first rigorous work characterizing some geometric properties of the ground state 
(global minimizer). 

We conclude this section on the interesting connection with gradient-flow dynamics. It is 



convenient to examine either the H~ gradient flow of (1.1) or the modified Mullins-Sekerka 
free boundary problem of Nishiura-Ohnishi [22] which results from taking the gradient flow of 
the sharp-interface functional. In |14[ [13] the authors explore the dynamics of small spherical 
phases (particles). By constructing approximations based upon an Ansatz of spherical parti- 
cles similar to the classical Lifshitz-Slyozov- Wagner theory, one derives a finite dimensional 
dynamics for particle positions and radii. Here one finds a separation of time scales for the 
dynamics: Small particles both exchange material as in usual Ostwald ripening, and migrate 
because of an effectively repulsive nonlocal energetic term. Coarsening via mass diffusion 
only occurs while particle radii are small, and they eventually approach a finite equilibrium 
size. Migration, on the other hand, is responsible for producing self-organized patterns. For 
large systems, kinetic- type equations which describe the evolution of a probability density are 
constructed. A separation of time scales between particle growth and migration allows for 
a variational characterization of spatially inhomogeneous quasi-equilibrium states. Heuristi- 
cally this matches our findings of (a) a first order energy which is local and essentially driven 
by perimeter reduction, and (b) a Coulomb-like interaction energy, at the next level, respon- 
sible for placement and self organization of the pattern. Moreover, in [13] . one finds that both 
the particle position radii and centre ODE's have gradient-flow structures related to energies 
which can be directly linked to our first and second order limit functionals, respectively. 

The natural question is to what extent one can rigorously address the dynamics and the 
separation of coarsening and particle migration effects. Recently, Niethammer and Oshita 
[21] have given a rigorous derivation of the mean-field equations associated with the evolution 
of radii. Another approach (currently in progress) is via Sandier and Serfaty's connection 
between T-convergence and an appropriate (weak) convergence of the associated gradient 
flows [291 131] . Le [TB] has recently used this framework for the e — > problem, establishing 
convergence of the fT _1 -gradient flow of (1.1) to that of the modified Mullins-Sekerka free 



boundary problem of Nishiura and Ohnishi |22j . While this method gives a rather weak 
notion of convergence, it allows for much weaker assumptions on the initial data and generic 
structure of the evolving phases. 
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